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Abstract: We compare mj-2 with itict] both are kinematic variables designed to find 
relationships between masses of pair-produced new states with symmetric decay chains. 
We find that for massless visible particles tuct equals m,T2 m a particular limit. We 
identify advantages and disadvantages to the use of each variable. Tovey's paper on itict 
also introduced a powerful concept of extracting mass information from an analysis at 
intermediate stages of a symmetric decay chain. We suggest that rriT2 is a better tool 
for performing this analysis than mcT due to m^'s better properties under initial state 
radiation. 
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Dark matter's likely signature in a hadron collider will be missing transverse momen- 
tum. The stability of dark matter suggests a charge or conservation law that requires dark 
matter particles be produced in pairs at colliders. The kinematic variables m,T2 introduced 
by Lester and Summers [l]and mcT introduced by Tovey [2] aid in the task of determining 
the mass of new states that decay to dark matter particles at hadron colliders. Although 
my 2 has been used extensively (see [3, 4, 5, 6, 7] for a few examples), the variable ttict is 
new but shares many similarities and differences with rriT2- This note briefly defines mT2 
and mcT, explains when they give identical results, when they differ, and comments on 
benefits of each in their intended applications. 

Both variables assume a pair-produced new-particle state followed by each branch 
decaying symmetrically to visible states and dark-matter candidates which escape detection 
and appear as missing transverse momentum. Fig 1 is the simplest example on-which we 
can meaningfully compare the two kinematic quantities. The figure shows two partons 
colliding and producing some observed initial state radiation (ISR) with four momenta g 
and an on-shell, pair-produced new state Y. On each branch, Y decays to on-shell states 
X and v\ with masses mx and m Vl , and X then decays to on-shell states N and vi with 
masses and m V2 . The four- momenta of vi, V2 and N are respectively ati, «2 and p on 
one branch and fli, 02 and q in the other branch. The missing transverse momenta is 
given by the transverse part of p + q. 

First we describe rriT2- The variable rriT2 accepts three inputs: x (an assumed mass of 
the two particles carrying away missing transverse momenta), a and (5 (the visible momenta 
of each branch), and fx = (p + q)t (the missing transverse momenta). The variable rriT2 
is the minimum mass of the pair of parent particles compatible with the observed particles' 
four momenta and an assumed mass for particles carrying away the missing momenta. We 
can define rriT2 m terms of the transverse mass of each branch where we minimize the 
maximum of the two transverse masses over the unknown split between p and q of the 
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Figure 1: A simplest topology with which we compare tot2 and mcT- 
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overall missing transverse momenta: 



"4 2 (XiV,a,/3, Pt) = mki [max {m|(a,p), m^(f3, q)}] . (1) 

In this expression xn is the assumed mass of N, a and (3 are the four momenta of the 
visible particles in the two branches, the transverse mass is given by m^(a,p) = m\ + 

X% + 2(Et(p) Et(ch) —pr ■ cut) and the transverse energy Et(p) = \jp\ + Xyv ^ s determined 
from the transverse momentum of p and the assumed mass of the particle associated with 
momentum p. An analytic formula for the case with no transverse ISR can be found in the 
appendix of [4]. For each event, the quantity itit2(xn = mjv> «i + &2, Pi +P2, Pt) gives the 
smallest mass for the parent particle compatible with the event's kinematics. Under ideal 
assumptions, the mass of the parent particle Y is given by the end-point of the distribution 
of this rriT2 parameter over a large number of events like fig 1 . Because a priori we do not 
know ttitv, we need some other mechanism to determine 1 . We use x to distinguish 
assumed values of the masses {xy, Xx, Xn) from the true values for the masses (my, mx , 
mjv). Because of this dependence on the unknown mass, we should think of maxmT2 &s 
providing a relationship or constraint between the mass of Y and the mass of N. This 
forms a surface in the {xy, Xx, Xn) space on which the true mass will lie. We express this 
relationship as Xy(xn) 2 • 

Tovey [2] recently defined a new variable mcT which has many similarities to m,T2- 
The variable is defined as 

mlriaufa) = (E T ( ai ) + E T (Pi)f - («it - Pit) 2 - (2) 

Tovey's goal is to identify another constraint between masses in the decay chain. He 
observes that in the rest frame of Y the momentum of the back-to-back decay products X 
and v\ is given by 

(Um Y ,m x ,m Vl )) 2 = ^ " ^ + " — " — ( 3 ) 

where k* is the two-body excess momentum parameter (2BEMP) 3 . In the absence of 
transverse ISR (gx = 0) and if the visible particles are effectively massless (m Vl = 0), Tovey 
observes that max mcr(ai,/3i) is given by 2k*; this provides an equation of constraint 
between my and mx- Tovey observes that if we could do this analysis at various stages 
along the symmetric decay chain all the masses could be determined. 

The big advantage of m-cr is in its computational simplicity. Also, mcr is intended to 
only be calculated once per event instead of at a variety of choices of x- I n contrast, mT2 



lr rhe true mjv and my can be found in the case where Y undergoes a three-body (X is off-shell) through 
kinks in tot2 [5, 6, 8] or when combined with endpoints from other distribution (like maxm;;) [7]. 

2 In principle this surface would be considered a function of \y (xx, Xn), but niT2 makes no reference to 
the mass of X and the resulting constraints are therefore independent of any assumed value of the mass of 
X. 

3 Tovey refers to this as the 2-body mass parameter Mi. We feel calling this a mass is a bit misleading 
so we are suggesting 2BEMP. 
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is a more computationally intensive parameter to compute; but this is aided by the use of 
a shared repository of community tested C++ libraries found at [9]. 

How are these two variables similar? Both mcT an d TriT2> m the absences of ISR, 
are invariant under back-to-back boosts of the parent particles momenta [6]. The variable 
mcT equals rriT2{x = 0) m the special case where x = and when the visible particles are 
massless (a 2 = f3f = 0) and the there is no transverse ISR (gx = 0) 

mcr(«i,/9i) = m T2 {x = 0,ai,/?i,^ T = (p + q + a 2 + (3 2 )t) if a? = fif = 0. (4) 
= 2(a 1T ■ Pxt + I«it| \Pit\)- ( 5 ) 

The mcT side of the equation is straight forward. The rriT2 side of the expression can be 
derived analytically using the formula for m,T2 given in [4] ; we also outline a short proof in 
the appendix. Eq(4) uses a rriT2 m a unconventional way; we group the observed momenta 
of the second decay products into the missing transverse momenta. In this limit, both 
share an endpoint of 2k* = {m\- — m 2 x )/mY- To the best of our knowledge, this endpoint 
was first pointed out by Cho et.al. [5] 4 . We find it surprising that a physical relationship 
between the masses follows from rriT2 evaluated at a non physical x- I n the presence of 
ISR, eq(4) is no longer an equality. Furthermore in the presence of the ISR, the end point 
of the distribution given by either side of eq(4) exceeds 2/c*. In both cases, we will need to 
solve a combinatoric problem of matching visible particles to their decay order and branch 
of the event which is beyond the scope of this paper. 

In the case where the visible particle v% is massive, the two parameters give different 
end-points 

m\-m\ m 2 vi 

maxmcr = 1 (6) 

my my 

maxm T2 (x = 0, fV = (p + Q + «2 + Pi)t) = J m\ + 2(&| + K^Jk 2 + mgj (7) 

where A;* is given by eq(3). Unfortunately, there is no new information about the masses 
in these two endpoints. If we solve eq(6) for mx and substitute this into eq(7) and (3), all 
dependence on my is eliminated. 

Tovey's idea of analyzing the different steps in a symmetric decay chain to extract the 
masses is powerful. Up until now, we have been analyzing both variables in terms of the 
first decay products of Y. This restriction is because txict requires no transverse ISR to give 
a meaningful endpoint. If we were to try and use a 2 and j3 2 to find a relationship between 
mx and m^, then we would need to consider the transverse ISR to be {g + a\ + (3i)t which 
is unlikely to be zero. 

We suggest is a better variable with which to implement Tovey's idea of analyzing 
the different steps in a symmetric decay chain because of its ISR properties. With and 
without ISR, rriT 2 s endpoint gives the correct mass of the parent particle when we assume 

4 The endpoint given by Cho et.al. is violated for non-zero ISR at \n < wiat and \n > iun- 
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Figure 2: Shows constraints from max mj-i used with different combinations as de- 
scribed in eqs(8,9,10) and the maxm^ described in eq(12). Intersection is at the true mass 
(97GeV, 144 GeV, 181 GeV) shown by sphere. Events include ISR but otherwise ideal conditions: 
no background, resolution, or combinatoric error. 



the correct value of the missing-energy-particle's mass . For this reason, maxmr 2 gives a 
meaningful relationship between masses (my,mjf,mjv) for all three symmetric pairings of 
the visible particles across the two branches. A relationship between my and mx is given 
by 

Xy{xx) = maxm T2 (xx,ai,/3i,-^T = (p + Q + «2 + 02}t)- (8) 
A relationship between mx and mjy can be found by computing 

Xx (xjv) = maxm T2 (xJV, "2, 02, fV = ip + ?)t) (9) 

where we have grouped a,\ + 0\ with the g as a part of the ISR. A relationship between 
my and mjy can be found by using rriT2 i n the traditional manner giving 

Xy{xn) = maxm T2 (xN,ai + a 2 ,0i + 02, fr = (p + q)t)- (10) 

Lastly, we can form a distribution from the invariant mass of the visible particles on each 
branch m\ 2 = (a± + a 2 ) 2 or m\ 2 = (Pi + 02) 2 ■ The endpoint of this distribution gives a 
relationship between my, mx, and m^ given by 

maxm 2 2 = ^ Y -ml){ml-m%) _ 
m x 

In principle we could plot the maxmT2(xx,ai,/3i,l l T = (a.2 + fiz + p + q)r) vs \x as a function of 
transverse ISR and the value of \x at which the end point is constant would give the correct value of mx', 
at which point the distributions end point would give the correct my- In practice we probably will not 
have enough statistics of ISR events. 
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Solving this expression for my gives the relationship 



Xy(xn,Xx) 



Xy((maxmf 2 ) + Xx ~ X%) 
x\ ~ Xn 



(12) 



Fig 2 shows the constraints from eqs(8,9,10,12) in an ideal simulation using (my = 181 GeV, 
rax = 144 GeV, mx = 97 GeV), 1000 events, and massless visible particles, and ISR added 
with an exponential distribution with a mean of 50 GeV. These four surfaces in principle 
intersect at a single point (my, mx, mpi) given by the sphere in the figure 2. Unfortunately, 
all these surfaces intersect the correct masses at a shallow angles so we have a sizable 
uncertainty along the direction of the sum of the masses and a tight constraints in the 
perpendicular directions. In other words, the mass differences are well-determined but 
not the mass scale. From here one could use a shape fitting technique like that described 
in [7] to find a constraint on the sum of the masses. Tovey's suggestion for extracting 
information from these intermediate stages of a symmetric cascade chain clearly provides 
more constraints to isolate the true mass than one would find from only using the one 
constraint of eq(10) as described in [5]. However, Tovey's suggestion is more feasible using 
the rriTi rather than ijict because the constraint surfaces derived from rriT2 intersect the 
true masses even with ISR. 

In summary, we have compared and contrasted mcr with mn- The variable mcT 
is a special case of rriT2 given by eq(4) when ISR can be neglected and when the visible 
particles are massless. In this case, the end-point of this distribution gives 2k*, twice the 
two-body excess momentum parameter (2BEMP). If m vi ^ 0, the two distributions have 
different endpoints but no new information about the masses. In the presence of ISR 
the two functions are not equal; both have endpoints that exceed 2k*. Because of it's 
better properties in the presence of ISR, rriT2 is a better variable for the task of extracting 
information from each step in the decay chain. Extracting this information requires solving 
combinatoric problems which are beyond the scope of this paper. 
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Appendix: Verifying in eq(5) 

We derived the side of eq(5) by following the analytic solution given by Barr and 
Lester in [4]. In this appendix, we outline how to verify that m,T2 is is indeed given by 



m T 2(x = 0,at,P,f> T 



cut — Pt) = 2(otT • Pt + \olt\ |/3x|) 



(13) 
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when a 2 = P 2 = and p 2 = q 2 = x 2 = and qt = 0. To do this we note that tot2 can 
also be defined as the minimum value of (a + p) 2 minimized over p and q subject to the 
conditions p 2 = q 2 = \ 2 (on-shell missing energy state), and (a + p) 2 = (/? + q) 2 (equal 
on-shell parent-particle state), and (a + P + p + q + g)x = (conservation of transverse 
momentum) [7]. 

The solution which gives eq(13) has pr = ~Pt and qx = — «t with the rapidity of p(q) 
equal to the rapidity of a ((3). We now verify that this solution satisfies all the constraints 
listed above. Transverse momentum conservation is satisfied trivially: (a + f3 + p + q)x = 
(a + (3 — a — (3)t = 0. The constraint to have the parent particles on-shell can be verified 
with 2|ar||pr| — ^Pt ■ c?r = 2|/?Tlkr| — ^Qt ■ Pt = 2|/3t||/9t| + 2<3t ■ Pt- 

Now all that remains is to show that the parent particle's mass is a minimum with 
respect to ways in which one splits up the missing transverse momentum between pr 
and qT while satisfying the above constraints. We take p and q to be a small deviation 
from the stated solution px = ~Pt + $T an d qT = -«t — <5t where St is the small 
deviation in the transverse plane. We keep p and q on shell at x = 0- The terms p , 
Pzi Qo, qz are maintained at their minimum by keeping the rapidity of p and q equal to 
a and (3. The condition that the parent particles are on-shell and equal is satisfied for 
a curve of values for St- The deviation tangent to this curve near \5t\ = is given by 
<5t(A) = A z x (aT\f3r\ + /3t|«t|) where x is a cross product, z denotes the beam direction, 
and we parameterized the magnitude by the scalar A. Finally, we can substitute p and q 
with the deviation 5r(A) back into the expression for the parent particle's mass (a + p) 2 
and verify that 2(ot • Pt + \&t\ \Pt\) at A = is indeed the minimum with respect to 
changes in A. 
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